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The free energy of inserting a protein into a membrane is determined by considering the variation 
in the spectrum of thermal fluctuations in response to the presence of a rigid inclusion. Both 
numerically and through a simple analytical approximation, we find that the primary effect of 
fluctuations is to reduce the effective surface tension, hampering the insertion at low surface tension. 
Our results, which should also be relevant for membrane pores, suggest (in contrast to classical 
nucleation theory) that a finite surface tension is necessary to facilitate the opening of a pore. 
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Bilayer membranes are self-assembled thin fluid sheets 
of amphiphilic molecules. They are characterized by 
small bending and large compression moduh, whose ef- 
fective values are influenced by thermal fluctuations . 
The softness of the bending modes permit large shape 
deformations which are important for the biological ac- 
tivities of some living cells (e.g., the red blood cell) 0- 
Biological membranes are typically highly heterogeneous: 
they usually consist of mixtures of different lipids and, 
in addition, contain a variety of different proteins which 
carry out diverse tasks such as anchoring the cytoskele- 
ton, opening ion channels, and cell signaling 

Membrane inclusions can modify the thermal fluctua- 
tions of the membrane by perturbing the local structure 
of the lipid matrix. It is well-known that the restrictions 
imposed on the thermal fluctuations of the membrane are 
the origin of attractive van der Waals-like forces between 
inclusions While these long-range interactions are 
typically very small, they are believed to play an impor- 
tant role in determining the phase behavior (e.g. aggre- 
gation) of such systems. Perturbing the spectrum of ther- 
mal fluctuations is also expected to contribute to the free 
energy associated with the insertion of proteins into lipid 
bilayers. This has an influence on the solubility of pro- 
teins and other membrane inclusions. Remarkably, this 
important entropic contribution to the insertion free en- 
ergy of a single protein has been ignored in previous cal- 
culations In this letter we study the free energy cost 
of inserting a rigid inclusion into a membrane, explicitly 
taking into account effects due to membrane fluctuations. 
These effects depend only on the inclusion's characteris- 
tic size. For transmembrane proteins, the magnitude of 
the fluctuation free energy can be as large as 25/cbT. As 
this contribution is comparable, but of opposite sign, to 
other free energy components, it strongly influences the 
thermodynamic stability of proteins. Our results should 
also be relevant for the fluctuation spectrum and nucle- 
ation energy of a membrane pore. At low tension, the 
fluctuation free energy acts as a barrier to the opening of 
a pore. 

We consider a bilayer membrane consisting of A'' lipids, 
that spans a planar circular frame of a total area Ap = 



in which a rigid inclusion of radius rg ^ Lp has 
been inserted. The Hclfrich energy (to quadratic order) 
for a nearly-flat membrane in the Monge gauge is given 
by i 



Hi = (jA.„ + - d 



(1) 



where cr is the surface tension, k the bending rigidity, and 
h the height of the membrane above the frame reference 
plane. The boundaries of integration in Eq.J^I include 
the outer (frame, r = Lp) boundary and the inner (in- 
clusion, r — ro) edge. The Laplacian in the Helfrich 
energy requires that we have two boundary conditions 
(BCs) for each boundary. On the inner boundary we fix 
the height of the membrane h{ro) = LL{(t>) and the con- 
tact slope dh{ro)/dr = H'{(f>), where (j> is the polar angle 
measured from the inclusion's axis of symmetry. On the 
outer boundary will have the natural BCs: h{Lp) — 
and V^h{Lp) — 0. The particular choice of outer BCs 
does not modify the free energy of the system in the 
thermodynamic limit. 

To gain insight into the contribution of thermal fluc- 
tuations to the insertion free energy we write the height 
function as h = ho + f where ho is the extremum of 
Hamiltonian i.e.. 



(2) 



subject to the BCs that /lo(ro) — H{(f), dh{ro)/dr ~ 
H'{(j)), ho{Lp) = 0, and V^ho{Lp) = 0. This implies 
/(''o) = and df{ro)/dr = on the inner boundary, and 
f{Lp) = and \7'^f{Lp) = on the outer boundary. The 
Helfrich energy can be written as 



ni (ho + /) = aAp 



a{Vhof + K{V^hoY (3) 



[aVho ■ V/ + KW^hoW^f] + ^[<J (V/)' + K (VV)' 



For the cross term (third term in Tii) we obtain, upon 
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integration by parts, 

/ dV [-crV^/io + kV^/io] / + / kV^/io (n • V) / + 

f (n • V) [^/lo - «V2/io] /, (4) 

where the last two integrals in the above equation are per- 
formed on the boundaries of the system, and n is a unit 
vector normal to the boundaries. The boundary terms in 
Eq.Q vanish since / = and fi ■ V/ = —df/dr = on 
the inner boundary, and / = and V^ft-o on the outer 
boundaries. The bulk term also vanishes by virtue of 
Eq.@. 

Without the cross term in Eq.Q, we are left with 
three terms: the projected area term aAp, the equilib- 
rium term depending on ho, and the fluctuation term 
depending on /. Thus, the energies associated with ho 
and / completely decouple and their contributions to the 
free energy are additive. Note that in our approach the 
equilibrium part of the free energy includes a contribu- 
tion from the height and tilt fluctuations of the inclu- 
sion. It is obtained by calculating the dependence of ho 
on the boundary values H{(j)) and H'{(j)), and perform- 
ing an appropriate thermal average over these quanti- 
ties. Other energetic components, such as hydrophobic- 
ity, translational entropy, electrostatics, should be added 
to the equilibrium term, and can be included in its defini- 
tion . The equilibrium term has been analyzed in many 
previous studies 5] . Its magnitude is protein specific and 
is usually in the range of —5 to —20kBT 8]. In contrast, 
the effect of membrane fluctuations on the insertion free 
energy has not yet been considered in the literature. We 
proceed to calculate the fluctuation part of the insertion 
free energy. Note that it is independent of the height 
and the contact angle of the inclusion (and their thermal 
fluctuations), which affect only the equilibrium part |^. 

Neglecting the equilibrium term, we are left with the 
projected area and the fluctuation terms. By integrat- 
ing the latter by parts twice, the remaining Hamiltonian 
takes the form 



n2{f) ^<yAp + ^ Jd^rf i-aW^ + kW^) f. 



(5) 



The boundary terms vanish in the above expression due 
to our choice of BCs: /(rg) = 0, df{ro)/dr = 0, f{Lp) — 
0, and V'^f{Lp) = 0. We proceed by expanding the func- 
tion / in a series of eigenfunctions fm nif) of the opera- 
tor C = -aV2 -t- kV*: /(r, 0) = /l„^,„/™,„(r)e™'^. 
The functions frn,n{r) can be written as the linear com- 
bination of the Bessel functions, Jm{r) and F,„(r), of the 
first and second kinds of order m, and the modified Bessel 
functions of the first and second kinds of order m, Km{r) 
and Im{r): 

+ CKrn{XT"r)+DIrniXT"r), 



equation Jm(AiLp) = 0. 
the limit that A^'^ro <C 



where the Xi {i = 1,2) are the positive solutions 
of (-l)^+V(Ar'")^ + f^iXm' = Mm,n, and ^„,„ is 
the eigenvalue corresponding to the function fm,n{r): 

Applying the BCs at tq and Lp, we derive the eigen- 
value equation 

Al [Im{hro)Kjn{X2Lp) - lm{hLp)K,n{X2ro)] 

K.(Airo)J,„(Aiip) - j;,(Airo)y,„(Airo)] 

= A2 [i^:„(A2ro)/„(A2ip) - /;,(A2ro)i^„(A2£p)] 

[Jm{Xiro)Yrn{XiLp) - Jm{XlLp)Ym{Xiro)] (6) 

(for brevity, we have omitted the superscript (m, n) from 
the notation of the A^ in the above equation). In contrast, 
for membranes without inclusions, we solve the simple 

It is interesting to note that, in 
m|, Eq.® reduces to the eigen- 
value equation in the absence of inclusions. This has 
the physically appealing interpretation that modes with 
characteristic lengths much larger than the inclusion ra- 
dius are hardly perturbed by its presence. In the opposite 
limit, A™'"ro ^ |to| (which also implies A™'"Lp ^ I'Tij), 
we can neglect terms proportional to /m(A™'"Lp) (which, 
otherwise, become exponentially large) and replace the 
remaining Bessel functions by their leading order asymp- 
totic expressions. This gives, for A"''" ^ ct/k, the 
simple equation tan [A'"'" (Lp — ro)] = 1, and the solu- 
tions A7'" « [|m|/2 + n+ (-1)" 7r/4] n/{Lp - To). The 
physical interpretation of this result is that the inclu- 
sion acts like a hard wall for modes with characteristic 
lengths much smaller than its radius. The effective linear 
size of the membrane for these modes is reduced from Lp 
to Lp — ro and the eigenvalues in this regime increase by 
roughly a factor of Lp/{Lp — tq). Thus, the dominant 
effect of the inclusion on the short wavelength modes is 
to lower the density of contributing modes in "A-space" 
[Note that A7'"+' - A^'" = ^/(Lp - ro)]. 

When the function /(r,<^) = E„,„ /im,n/m,n(r)e*""^ 
is substituted in Hamiltonian we find, due to the 
orthogonality the eigenfunctions 



rdr/„i,„i(r)/„2,n2(r)e'("i+"2)^ 



flO Srnl,-m2 Snl,n2, 



(7) 



that the modes decouple and that the Hamiltonian takes 
a quadratic form in the amplitudes |ft.m,n|. The normal- 
ization coefficient aq in Eq.(|7| is the projected area per 
amphiphilic molecule in the bilayer. Tracing over |/im,n| 
leads to the following expression for the Gibbs free energy 
associated with Hamiltonian Ti2 



G{a,Ap) = 
2 ^ 



(8) 



In. 



[a(A7^")2 + k(X^^^Y\ ApX, 
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where AdB is the thermal de-Broglie wavelength of the 
lipids. The Helmholtz free energy is given by F(A, ^p) = 



3 



G — a A, where the total membrane area A is related to 
the surface tension by 



(a) 



A — A„ 



1 



(9) 



Assuming that the membrane is incompressible and, 
therefore, that its total area is fixed, we can use Eq.Q to 
derive the following equation, relating the surface tension 
and the inclusion's radius 



1 



1 



= 0. 



(10) 

In the above equation A^^g^ are the corresponding solu- 
tions of the eigenvalue equation in the absence of the in- 
clusion (ro = 0): Jm(A™(p^Lp) = 0, and cto = cr(ro = 0). 
The solution to Ea. (|10|l has the form 



(To(l + (5), where 5 OijQ/Lpf 



(11) 



The projected area and fluctuation parts of the insertion 
free energy Ai^(ro) = F(ro) — i^(0) can now be calculated 
using Eqs.® and ifTUI) . We find that AF(ro) is given by 



AF(ro) 



(12) 



(0) 



)2 + «(A™-)4 
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Note that only ctq appears in the above expression, which 
is due to Ea. lfTTl) and the fact that we attempt to calcu- 
late Ai^(ro) only up to quadratic order in ro/Lp. For the 
same reason we can use ctq rather than a in the eigenvalue 
equation l^. The surface tension appears implicitly in 
this equation, through the relation A^ = A^ -I- a / k. In 
expression (|12(l we assume that the number of molecules 
forming the bilayer membrane does not change with the 
insertion of the protein. Consequently, the total number 
of modes which is directly proportional to the number 
of molecules in the bilayer is kept constant. In contrast, 
the projected area per molecule [which appears in Eq.Q] 
does depend on the radius of the inclusion, and this is the 
origin of the term (L^ — r^)/ L-^ appearing in the argu- 
ment of the logarithm in Ea. (|12|l . The first term on the 
right hand side (r.h.s) of Eg. 1)12(1 comes from the reduc- 
tion of the projected area. We will now show that, to 
a good approximation, the second term on the r.h.s. of 
Ea. (|12|l is quadratic in rg and, thus, can be interpreted 
as a thermal correction to the surface tension. 

In order to obtain an analytical result for the free en- 
ergy H12() . we make the approximation [based on our 
discussion of the asymptotic behavior of the eigenval- 
ues A™'", see the text after Eq.®] that eigenvalues such 



that Aj ' ro < 



(long wavelength) are not affected 



by the inclusion, whereas modes with Aj^ ' rp > Q!|to| 
(short wavelength) grow by a factor Lp/{Lp — rg). The 
numerical constant a is of the order of unity and its 
value, which may depend on the surface tension ctoj will 
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FIG. 1: The insertion free energy /S.F as a function of the 
inclusion's radius for k = IQksT and various values of (tq. 
The inset to graph (a): a log-log plot of the numerical results 
for (Jo = 0. The slope of the straight dotted line is 2. 



be determined later by an exact numerical evaluation of 
/S.F. We have verified numerically that this asymptotic 
form is indeed correct. We set n = 0, 1, ... , v'TVo, and, 
m = — -\/7Vo, . . . , \/Nq so that the total number of modes 
(degrees of freedom), 27Vo, is proportional to the num- 
ber of molecules forming the membrane, N . Along with 
these approximations, we evaluate the sum in equation 
(|12|l as an integral, giving us the simple result (correct 
up to quadratic order in tq) that A_F = — 7r((To — cr*)rQ, 
where 




(13) 



^ = k/cto, and £o = Lp/^/No is a microscopic length 
cutoff on the order of the characteristic size of a mem- 
brane molecule. We thus obtain the result that the fluc- 
tuations renormalize the surface tension. It is interesting 
to note that this renormalization tends to occur with the 
opposite sign as the bare surface tension (for £q < ^), 
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thus making it harder to insert an inclusion. Only for 
very stressed membrane ^ ^o) does a* become nega- 
tive. This is due to the reduction of the projected area 
that allows more thermal fluctuations. A more careful 
analysis of the long wavelength modes shows that these 
contribute only finite-size effects to the free energy which 
vanish in the limit of Lp 3> tq. 

We have numerically solved the eigenvalue equation ® 
and used the solutions to evaluate the sum in Ea. lfT^ . 
Numerical values of AF(ro) (for k = lOfc^T and various 
values of (Tq) are shown in Fign(a)-(b). They have been 
extracted by extrapolating the numerical results obtained 
for several values of 750 < Nq < 2000 to the thermody- 
namic limit A'o oo. In the inset to Fig^ (a), the 
results for ctq = are replotted on a logarithmic scale, 
showing that our prediction of a quadratic relation be- 
tween AF and tq is attained only for large inclusions 
with rg ^ 100^0 (the slope of the straight dotted line 
is 2). This is a typical size for colloidal particles [T^ . 
The value of the constant a appearing in Ea. lfT^ shows 
a slight dependence on the surface tension varying from 
1.59 for (To = to 1.72 for ^ = ^k/cto = 5?o/vr. The 
solid curves in Fign(a)-(b) depict our analytical expres- 
sion for AF, with a determined by fitting the results for 
large to Ea. ll3|l . From Fig^ (a) we conclude that, 
because of thermal fluctuations, there is a free energy 
penalty to embedding an inclusion in a weakly stretched 
membrane (small cto). For transmembrane proteins with 
typical radii of ro ^ 5£o, the energy cost is AF < 25fcBT, 
which is comparable to the equilibrium contribution but 
of opposite sign. This demonstrates the importance of 
the membrane fluctuations in determining the distribu- 
tion of transmembrane and free proteins. For larger in- 



clusions, the fluctuation free energy will dominate the 
equilibrium part. On the other hand, Fig^ (b) shows 
that inclusions greatly reduce the free energy of strongly 
stretched membranes (large cfq). The primary reason 
that the free energy is lowered in this regime is the re- 
duction of the projected area. These results should also 
be relevant for the question of nucleation of a membrane 
pore which, albeit more complicated, can be studied by 
similar approach 0. They suggest that there exists a 
(finite!) critical value of the surface tension below which 
pores cannot open and above which they grow without 
bounds. Classical nucleation theory, which ignores fluc- 
tuations effects, predicts that the critical surface tension 
is zero Q. 

In summary, we have computed the free energy of in- 
serting an inclusion into a membrane. We explicitly cal- 
culated the contribution of membrane fluctuations. The 
primary effect of these fluctuations is to reduce the effec- 
tive value of the surface tension. At low surface tension 
it provides a positive component to the free energy of 
an embedded inclusion, thereby impeding the insertion 
of transmembrane proteins. The sensitivity of the free 
energy to variations of the surface tension suggests that, 
by controlling the membrane surface tension appropri- 
ately, one may control the thermodynamic stability of 
embedded proteins and, thus, the equilibrium distribu- 
tion between proteins inserted in the membrane and in 
solution. 
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